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ABSTRACT 

We perform two-dimensional hydrodynamical simulations to qua ntitatively explore the torque bal¬ 
ance criterion for gap-opening (as formulated bv iGrida et aUl^OOdH in a variety of disks when consid¬ 
ering a migrating planet. We find that even when the criterion is satisfied, there are instances when 
planets still do not open gaps. We stress that gap-opening is not only dependent on whether a planet 
has the ability to open a gap, but whether it can do so quickly enough. This can be expressed as 
an additional condition on the gap-opening timescale, tgap, versus the crossing time, Grossj be. the 
time it takes the planet to cross the region which it is carving out. While this point has been briefly 
made in the previous literature, our results quantify it for a range of protoplanetary disk properties 
and planetary masses, demonstrating how crucial it is for gap-opening. This additional condition has 
important implications for the survival of planets formed by core accretion in low mass disks as well 
as giant planets or brown dwarfs formed by gravitational instability in massive disks. It is particu¬ 
larly important for planets with intermediate masses susceptible to Type Ill-like migration. For some 
observed transition disks or disks with gaps, we expect that estimates on the potential planet masses 
based on the torque balance gap-opening criterion alone may not be sufficient. With consideration of 
this additional timescale criterion theoretical studies may find a reduced planet survivability or that 
planets may migrate further inwards before opening a gap. 

Subject headings: hydrodynamics — methods: numerical — planet disk interactions — planets and 
satellites: dynamical evolution and stability — protoplanetary disks 


1. INTRODUGTION 

To predict the evolution of a planet’s orbital separa¬ 
tion in a gas disk it is crucial to know the conditions 
leading to different types of migration. Moreover, for a 
planet to survive the lifetime of the disk (rs 1 — 10 Myr; 
iHaisch et al.l I200ll : lMamaie^l2009fl without falling into 
the central star the radial motion of the planet needs 
to be slow enough. This can be achieved, for exam¬ 
ple, through the formation of a gap around the planet, 
which through a deficit of disk material in the planet’s 
co-orbital region decreases the angular momentum ex¬ 
change between the planet and the disk. 

1.1. Gaps in observations 

Annular regions with locally lower surface mass den¬ 
sities have been inferred based on the analysis of 
spectral energy distribut ions (SEDs) for decades (e.g. 
IMarsh &: Mahone^ 119921) . However, these interpreta¬ 
tions of the data are ambiguous given uncertainties in 
the opacity and temperature distributions. In recent 
years, inner holes (e.g. CoKu Tau 4) and gaps (e.g. 
GM Auriga) have been confirmed by detailed high signal- 
to-noise SEDs fro m the Spitzer Space Telescope (e.g. 
ICalvet et al.ll200^ . Inner holes can be understood ei¬ 
ther as a natural consequence of disk evolution due to 
grain growth (and accompanying changes in opacity) or 
photoevaporation ([Alexander et al.ll2013D or as an indi- 
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cation for the presence of planets. As the midplane tem¬ 
perature falls gradually in circumstellar disks, gaps are 
hard to identify from SEDs unless they are very large. 
Gaps are more often identified by direct imaging, either 
in scattered light or thermal emission. One caveat to 
these interpretations is that the appearance of a gap in 
scattered light could also be explained by shadows in the 
outer disk due to structures in the inner disk. In most 
cases of resolved gaps, the innermost disk is unresolved, 
but inferred from SEDs and confirmed by on-going ac¬ 
cretion onto the star. The inner edge of t he outer disk 
is oft en > 20 AU and can be resolved (e.g. iGarufi et al.l 
l2013f) . Thus there is a growing list of objects with re¬ 
solved outer disks and strong evidence for an inner disk 
with gaps ^ 10 AU or more. There are far fewer ob¬ 
jects with direct observational evidence for gaps where 
the outer edge of the inner disk and the inner edge of the 
outer d isk are both resolv ed in multiple bands (e.g. HD 


100546 iWalsh et al 


169142 iQuanz et al 


20M & iPineda et al.l 120141 and HD 


201 , 1 ). 


1.2. Theory of gap formation 

In the past several criteria have been proposed that 
quantify the planet and disk conditions leading to gap 
for mation. A pressure stab ility condition was proposed 
byE in fc PapaloizorJ (|1993D and is given by 

H < Rh ^ Rpiq/Sy/^ (1) 

stating that the pressure scale height H in the disk at the 
planet’s location, Rp, should be smaller than the planet’s 
Hill radius, Rh, for a gap to form. As Rh depends on the 
planet to primary star mass ratio, q, a sufficiently large 
planet mass is needed for gravitationally induced density 
wakes to result in shocks repelling the disk material from 
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the planet. This condition is also known as the ’’strong 
shock limit” or the thermal cr i terion for gap opening. 
Furthermore iLin fc Papaloizoul l)1993li stated that for a 
gap to form the disk viscosity should be small enough 
so that the disk’s viscous diffusion does not completely 
negate the mechanism of gap clearing. This condition on 
the viscosity can be written as 


Mp> 


dOr-pM* 


( 2 ) 


where M* is the primary mass, rip, Vp and Mp are the 
angular Keplerian velocity and the kinematic viscosity at 
the planet’s location and the planet’s mass, respectively. 
T he most general cri terion for gap-opening is provided 
by ICrida et al.l (|2006f) . 


3_ff 50^^^ 
4 Rh q^pRl - 


( 3 ) 


It is a semi-analytical criterion based on the balance be¬ 
tween pressure, gravitational and viscous torques for a 
planet on a fixed circular or bit. Hence we refe r to it as 
the torque balance criterion. ICrida et al.l (j200()[) define a 
gap when the mass density drops to 10% of the unper¬ 
turbed density at the planet’s location. 

In addition to the effects of the pressure and vis¬ 
cosity which act against gap formation, there i s also 
the time aspect to consider. iHourigan fc Ward! ()1984[ 1 
stated that the time required for gap clearing, tgap^ 
must be smaller than the planet’s migration timescale 
across its own horseshoe region, tcross, i-e. tcross ^ tgap- 
IWard fc HouriganI (|1989tl refer to this criterion as the 
’’inertial mass limit” given by 


®imit 



( 4 ) 


as this inequality imposes a lower limit for the plane¬ 
tary mass required for gap-opening. Ep oc R~^ repre¬ 
sents the unperturbed surface mass density that would 
be found at the planet’s location. It should be noted 
that these timescale considerations were calculated for 
low planetary masses in low-viscosity disks. Hence 
the problem could be solved analytically using linear 
approximati ons for the tidal t orque between the disk 
and planet. iLin fc Papaloizoul (1198611 studied the same 
timescale problem numerically and found the less strict 
limit tcross ^ = iH/R)tgap by arguing that it is not 

necessary to open a full gap during the crossing time but 
it suffices to significantly perturb the surface density to 
considerably alter the planet’s migration regime. 

Comparing the pressu re stability condit i on on the 
planet’s mass (eq. [J) and IHourigan fc WardI (Il984ll ’s in¬ 
ertial mass limit (eq. |4]) one finds that for a disk with 
profile fc = 3/2 the pressure stability condition is stricter 
if 


both criteria use linear analysis for the planet’s interac¬ 
tion with the disk which may break down when a planet 
grows in mass. 

We are curious that the inertial limit is dependent on 
the disk mass and surface mass density profile (eq. S]), 
but both parameters are missin g in the torque bal ance 
criterion (eq. [3]). In addition ICrida et al.l (|2006ll ne¬ 


glected planet migration. However, since these parame¬ 
ters directly affect the angular momentum exchange be¬ 
tween a planet and a disk, and hence planet migration, 
they may affect the planet’s ability to open a gap. 

In the last few years the pressure stability criterion 
(eq. [T|) or the torque balance criterion (eq. [3|) have 
been widely used to predict or explain numerical studies 
of planet formation/synthesis using various disk mod¬ 
els with some inves tigations of high-mass companions in 


massive disks (e .g. |AliberLel_aL [2004 [Ida^^J^ 


Zhu et all [Ml : 


Galvagni fc Maver 


Forgan fc Ricel 12014 IVorobvov 
1201411 . 


2004; 


2013; 


1.3. Outline of this paper 

As both the torque balance criterion and inertial mass 
limit are based on linear analyses of planet-disk interac¬ 
tions we are hence intrigued by the following questions: 

• Is the timescale condition still implicitly covered by 
the torque criterion in migration scenarios where 
non-linear torques prevail? 

• Can t he torque balance criterion of ICrida et al.l 
(1200611 be safely applied to disk models that do not 
necessarily offer the conditions usually assumed for 
a ’’pure” Type I migration? 

• Is gap-opening affected by the disk mass? 

We note that when considering gap-opening in disks 
(e.g. when interpreting observations) the criterion based 
on the balance of torques is always used but that based 
on the timescale is often not considered. In this work we 
point out that one needs to consider not just whether a 
planet can open a gap (based on the balance of torques) 
but also whether it can do it quickly enough in a vari¬ 
ety of disks expected to form companions by either core 
accretion or gravitational instability. 

As a broad distribution of disk and companion masses 
exists we choose two regimes for our studies: massive 
self-gravitating (SC) disks with migrating companions 
of the (sub-)brown dwarf regime and disks based on the 
Minimun Mass Solar Nebula (MMSN) model where we 
introduce Jupiter-mass planets. We let these companions 
migrate freely and compare our findings to the prediction 
of the torque balance criterion given by eq. [31 

In Section Fl.dl we outline our notation. We describe our 
numerical method and simulations performed in Sections 
[2] and [31 respectively. We then present our results in 
Section [H Finally, in Section [5] we discuss our results 
in the context of both theories and observations, and 
present our conclusions in Section [HI 


M* > 0.77ri?2Ep. 


( 5 ) 


For typical circumstellar disks this inequality is sat¬ 
isfied and thus planets satisfying the pressure crite¬ 
rion and consequently the stricter torque balance cri- 
terion o f ICrida et al.f (l2006f) will automatically satisfy 
IHourigan fc War*? ! 19841) ’s inertial mass limit. However, 


1.4. Notation 

We refer to our migrating objects as ’’companions” re¬ 
gardless of whether they would commonly fit into the 
brown dwarf or giant planet mass regime. 

The companion to primary mass ratio is dehned as 
q = Me/NR and the companion’s Hill radius is described 
by Rh ~ Rc{q/3y^^. The subscript ”c” indicates that 
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the parameter is evaluated for the companion or at the 
companion’s position. H = Cg/fl is the pressure scale 
height in the disk with the Keplerian an¬ 

gular velocity for a circular orbit and h = H/R is the disk 
aspect ratio. Finally, we write the physical disk viscosity 
V according to the a-disk model ly = ac.,H, with the 
dime nsionless stress parameter a (IShakura fc SunvaevI 

[Ml). 


2 . NUMERICAL METHOD 

The simulations are carried out using the polar 2D 
grid-based hydrody namics c ode, FARGO-ADSql, origi¬ 
nally developed by iMassetl (j200flf l and later extended 
to impleme nt a self-gravity solver base d on fast Fourier 
transforms (iBaruteau &: Mas^l2008bll . The disk’s grid 
is described by polar coordinates (i?, 0) with the star at 
the origin. An indirect term is included, which accounts 
for the frame acceleration due to the gravitational in¬ 
teraction of the bodies. The grid has open boundaries. 
The companions are represented by point masses and 
self-gravity is included in all of our simulations. 

2.1. Equations of state 

Our disk models include two distinctive regimes. For 
gravitationally unstable disks in which planets are ex¬ 
pected to form by gravitational instability (see section 
\Tni we use an adiabatic equation of state together with 
the energy equation 

^ + ^iult) = -p^ -It + Q+-Q., ( 6 ) 

where u is the thermal energy density, it is the flow 
velocity and p is the ve rtically integrated pressure 
([Baruteau fc: Mas^l2008ail . The first term on the right 
hand side describes the compressional heating due to 
physical processes and Q+ expresses the shock heating 
due to artificial bulk viscosity. Q- = u/tcooi is the cool¬ 
ing term with Uooi being the cooling timescale. The 
equation of state is written by 

P = (7 - l)u = YsRspecT. (7) 

where Rgpec = Runiv/ a) is the specific gas con¬ 
stant {Runiv is the universal gas constant, p = 2.4 
is the mean molecular weight of the gas particles and 
nip is the proton mass), Na is the Avogadro Number 
and T is the temperature. The sound speed, Cs, is 
given by cl = {dp/dH), which for the adiabatic model 
(p = const ■ T,'^) leads to Cg^adiab = yj 7P/U. The adia¬ 
batic index is 7 = 5/3. _ 

For MMSN-like disks (see section 13.2|) we use an 
isothermal equation of state without solving an energy 
equation. In these isothermal runs the viscosity is added 
in the momentum equation only. Here, the sound speed 
is given by Cs^iso = Cs^adiab!y/l- 

2.2. Gravitational instability 

The gravitationally unstable disks are set up so that 
the self-gravitating structure is present but the disks do 
not fragment (so that the interactions between just a 

® See fargo.in2p3.fr 


single companion and its disk can be considered). The 
fragmentation of d isks requires t wo conditions to be sat¬ 
isfied. Firstly the (|Toomrelll964D criterion should be ful¬ 
filled, i.e. 


Q 


CgO 

ttGS 


< 1 , 


( 8 ) 


where Q is the dimensionless Toomre parameter. The 
second condition for fra gmentation req uires the disk to 
cool on a fast timescale l)Gammiil200Ill . Expressing the 
cooling timescale as a dimensionless cooling parameter 
P = tcooi^ the critical value below which the disk frag¬ 
ments was estimated to be Peru ~ 3 and 7 for a ra¬ 
tio of specific heats of 2 and 5/3, r espectively (iGammid 
I 2 OOH iRice et ;iril2005h . Recently, iMeru fc Batel ( ^llH 
showed that the numerical simulations leading to the lat¬ 
ter result were limited by the resolution. Upon increasing 
the resolution they found that Peru > 20 and may even 
be as large as 30 as a larger resolution involves more 
heating from the artificial viscosity emplo yed in numeri¬ 
cal si mulations to model shocks correctly (|Meru fc Bate! 
2 OI 2 II. While this is still an active area of research (e.g. 
Rogers fc Wadsle'^l20I2llPaardekooDeiIl20ll. IRice et al.l 

2014 iPaardekooper et al.ll201lL iLodato fc Clarkd 1201 III 
we set the cooling rate per unit area to be Q- = u/tcooi = 
{uQ)/P, where P = 30, to ensure no fragmentation oc¬ 
curs, and the Toomre parameter stabilizes in these equa¬ 
tions to Q ~ 2 throughout the disk. 

If we assume that the transport of angular momentum 
due to gravitoturbulence occurs localljo we can connect 
the stress para meter a with the cooling parameter P by 
(|Gammiell2001h 


4 1 _ 4 1 

97 ( 7 - 1 ) 12 ^ 00 / 97 ( 7-l)/3’ 


This indicates that a more rapid cooling increases the 
viscous stress a in the disk as the self-gravity is less effi¬ 
ciently countered by heating. 


2.2.1. Mechanical equilibrium 

A SG disk is self-regulating in the sense that inter¬ 
nal cooling and shock heating due to gravitoturbulence 
leads to a Toomre Q parameter near unity. The disk 
material rearranges naturally until a mechanical equi¬ 
librium with dYj/dt = 0 is reached. Using the viscous 
evolution equation for Keplerian disks in circular mo¬ 
tion from |Lynde n-Bell fc Pringlel (jl974[ l this condition 
becomes (jPringlelll981f l: 








= 0 . 


( 10 ) 


In our set-up the SG disks rearrange to the state E oc 
R~^/^ and T = const, which solves the above equation 
when inserting n = aCgH and using the fact that the self¬ 
gravity mechanism leads to a uniform Toomre parameter, 
Q = const. 


^ This does not have to b e gen erally true as was po inted out by 
IBalbus Hi Papaloizoul II1999I 1 and ICossins et al.l II2009I 1. 
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2.3. Softening length 

2D polar disk models do not strictly simulate the phys¬ 
ical impact of the vertical dimension of a realistic disk. 
For this reason a softening length is introduced which 
mimics the dampening of the companion’s potential in a 
disk with finite thickness given by 


GMc 

(d2+£2)1/2’ 


( 11 ) 


where $c is the companion’ s potential and d is t he dis¬ 
tance from the companion (iKlev &: NelsonI [201211 . Fur¬ 
thermore the utilization of a softening length of the order 
of a significant fraction of the disk scale height avoids 
mesh singularities in numerical grid-based simulations. 
The self-gravity of the d isk is treated analo gously with 
the softening length esG- iMiiller et al.l (|2012ll extensively 
studied the effects of various softening lengths. In their 
Figure 4 they show which value of the softening length 
should be considered to obtain the best approximation 
of a realistic physical force at a given distance from the 
companion. At very close ranges the suggested value 
varies significantly with the distance from the compan¬ 
ion. However, beyond a distance ~ Q.7H from the com¬ 
panion the value of Cc ~ 0.7iJ remains below an error of 
less than 10% from the optimum value. Following this 
reasoning we choose Cc = 0.7iJ and analogously, based 
on their Figures 13 & 15 we choose esc = 0.8i7. 


3. SIMULATIONS 

In this section we describe the disk configurations used 
in our simulations. First, we focus on massive gravita¬ 
tionally unstable disks. Then we switch to the low-mass 
end of the common circumstellar disk distribution by set¬ 
ting up disks based on the MMSN model to examine 
whether our results hold for a variety of disks and com¬ 
panion masses. 


3.1. Massive SG disk models 


We let the self-gravity develop over 30 orbits (orbital 
timescale at 100 AU) in the reference SG disk and over 
90 orbits in the lighter SG dislQ so that the disk reaches a 
quasi-steady state without any companions. During this 
time the disk parameter profiles rearrange as described 
in subsection 12 .2.11 and the resulting parameters can be 
seen in Tabled Note that due to gravitoturbulence some 
of the parameters deviate somewhat from their values 
shown in Table [5] as they suffer from stochastic fluctu¬ 
ations. The value of Q at a certain radial separation is 
estimated by 


Q 


(Cs) 


azimuth 


(H) 


azimuth 


7rG(E) 


azimuth 


( 12 ) 


where (■■•)azimuth denotes the azimuthal average of the 
physical quantity. 

Note that for our numerical studies we use simplifying 
models of circumstellar disks and so their parameters are 
not necessarily of the form found through observations, 
but are chosen so that the results can be understood 
under controlled conditions. 


3.2. MMSN-like disk models 

For the simulations at the lower end of the disk 
mass distribut i on w e set-up disks similar to a MMSN 
disk (iHavashil Il981h and the disk model utilized by 
iGrida et al.l ( 200fih . i.e. with typical T Tauri disk proper¬ 
ties in which planets would be expected to form by core 
accretion. The disks extend between 2.5 - 30 AU with a 
surface mass density profile E oc ^nd a constant 

temperature across the whole disk. 

We model three disks with different masses all satisfy¬ 
ing the above conditions as follows: 

1. a disk with Mdisk = 0.02Mq, which we call herein 
the ’’lighter MMSN disk”. 

2 . a disk with Mdisk = O.O5M0, wich we call herein 
the ” intermediate MMSN disk”. 


We use two different disk set-ups for our massive SG 
disk simulations: 


3. a disk with Mdisk = O.OSMq, which we call herein 
the ’’heavier MMSN disk”. 


1. The initial disk condition by iBaruteau et al.l 
(|201ir ). We refer to this as the ’’reference” con¬ 
figuration. 

2. The same as (1) but with half the initial disk mass, 
herein after called the ” lighter SG disk”. 

The initial conditions for the two SG disk models are 
described in Table [1] Most of the parameters are defined 
directly in the code. Others are evaluated indirectly, e.g. 
the temperature T, which is found by inserting the ex¬ 
pression for the sound speed into eq. [71 solving for T 
and using the relation = hRVL. Further we estimate 
the a-parameter using eq. [9] with /? = 30. In general a is 
measured as the sum of the gravitational and Reynolds 
str esses in the disk . As w e employ a similar disk set-up 
as IBaruteau et al.l (|20lH) , we adopt their findings (see 
their Fig. 2) suggesting that the real value of a might 
deviate from theory (equation jS]) by up to factor of two. 
The implication of such a deviation is discussed in section 


The detailed initial parameters for the disks are found 
in Table |3l 

3.3. Predictions for gap-opening 

Fig. [Tlsh ows the prediction of the semi-analytical crite¬ 
rion (eq. |3l IGrida et'mi200fi[l in each of our disk set-ups 
in a companion mass versus orbital radius diagram. The 
solid red line represents the criterion’s limit, i.e. eq. [3] 
with an equality, calculated using the values in Tables 
[2] and [3] for the SG and MMSN-like disks, respectively. 
Gompanions situated above this line are expected to open 
a gap. Note that this criterion is not dependent on the 
disk mass so the graph for the MMSN-like disks is appli¬ 
cable to all the disk masses considered since the thermal 
and viscous properties remain the same. In the SG disks 
the thermal properties vary with disk mass, hence the 
need for separate predictions for different disk masses. 

^ Due to the lower disk mass the self-gravity mechanism develops 
slower and the mechanical equilibrium is achieved at later times in 
the lighter SG disk compared to the reference SG disk. 
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TABLE 1 

Initial parameters for the SG disk models. 


disk model 

reference SG disk 

lighter SG disk 

parameters common to both disks 

grid size 

number of grid cells 
softening lengths 

primary mass 

o-stress parameter 

aspect ratio 

temperature 

20 - 

1536 (azimuthal) 
Ec = 0.7H, 

M* = 

a S3 

h = 

T = 15 K 

250 AU 

X 516 (radial,[log]) 

= 0.8J7 

= 1M0 

0.013 

= 0.1 

( ^ i~^ 

V lOOAU ) 

disk specific parameters 

disk mass 

surface mass density 

Toomre Q-parameter 

A^disk = O.4M0 
S = 25 g cm“^ 

Q = 1.7 

f R ^ 

-2 

A^disk = O.2M0 

S = 12.5 g cm~^ 

Q = 3.3 

( R ^ 

-2 

lOOAU ) 

y lOOAU ^ 


TABLE 2 

Parameters for the SG disk models at the time of companion introduction, when the self-gravity has already developed. 


disk model 

reference SG disk 

lighter SG disk 

parameters at companion introduction 


aspect ratio 
temperature 
disk mass 

surface mass density 
Toomre Q-parameter 


i 0.11 


\ lOOAU / 


T ss 17 K 
JW'disk = 0.29 Mq 
E ss 18.3 g cm“^ 
Q Ri 1.6 


-3/2 


,1, ~ 0.06 jQOAu^ 
T S3 6.2 K 
.^disk = O.15M0 
S S3 9.6 g cm“^ 

Q Si 2.0 


-3/2 


TABLE 3 

Parameters for the MMSN-like disk models at the time of companion introduction. 


disk model 

lighter MMSN-like disk 

intermediate MMSN-like disk 

heavier MMSN-like disk 

common parameters 

grid size 

number of grid cells 
softening length 

primary mass 

o-stress parameter 

aspect ratio 

temperature 

2.5 - 30 AU 

1536 (azimuthal) x 508 (radial,[log]) 

€c = 0.777, esG = 0.877 

M* = 1M0 

a = 0.0013 

T = 129 K 

disk specific parameters 

disk mass 

surface mass density 

Toomre Q-parameter 

JW'disk = O.O2M0 

S = 338 gcm-2(5fp) 

Q = 18 

A7disk = O.O5M0 

S = 845 gcm-2 

Q = 7.0 

A^disk = O.O8M0 

S = 1352 g cm-2 

Q = 4.4 


We allow a single companion in each of our disk set¬ 
ups to migrate freely and compare its ability to open 
gaps with the pred iction of the torque balance criterion of 
iCrida et al.l (|2006ll in equation ((S]). To be consistent with 
their approach we regard a gap as a one order of mag¬ 
nitude decrease in the surface mass density in the com¬ 
panion’s horseshoe region compared to the unperturbed 
density. The crosses in Fig. [1] show the radial loca¬ 
tion where each companion is introduced and the dashed 
horizontal lines indicate the inward migration tracks. All 


companions are expected to start gap clearing eventually 
as they cross the criterion’s line because of their inward 
migration. 

To test the criterion using our reference SG disk model 
we choose 15, 30, 45 and 80 Mjup for the companion 
masses. While these values a re higher th a n the initial 
clump mass estimates of e.g. IBolev et'HI (|2010ll given 
our disk parameters, clumps may accrete gas and hence 
companions formed by gravitational instability could be 
much more massive during their migration than at the 
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Gap prediction, reference SG disk 



Gap prediction, lighter SG disk 
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Gap prediction, MMSN-like disks 



Fig. 1. — Companion’s mass versus orbital radius indicating the torque balance criterion of ICrida et al.l II2006I) (solid red line). Companions 
are expected to open a gap when situated above the limit. The crosses indicate where in the disk the companions are introduced in our 
simulations while the dashed lines indicate the inward migration tracks. Left: Prediction for the reference SG disk. Middle: Prediction 
for the lighter SG disk. Right: Prediction for the three MMSN-like disks. All simulations are expected to show gap-opening at some stage 
during an inwards migration. 


time they form (e.g. IZhu et ani2012tl . 

Furthermore, we want to stress that the purpose of this 
study is a numerical experiment to test the gap-opening 
criterion. We therefore choose companion masses that 
are both above and below the critical mass limit given 
by the red line in Fig. [U left panel. To highlight the 
significance of the disk mass on the simulation results 
we utilize the same companion masses in our lighter SG 
disk. In the MMSN-like disks we employ 0.5, 1 and 2 
Mjup mass companions with the same intention. 

3.4. Introduction of the companions 

First, the companions are inserted into the SG disks at 
an orbital separation of 140 AU. A companion is intro¬ 
duced with a linear mass accumulation over one orbital 
period at 100 AU for the massive SG disk simulations (in 
Sect. 13.11) and one orbital period at 5 AU for the MMSN- 
like disk simulations (see Sect. 13.21) . Since the compan¬ 
ions are not fixed to a particular radial separation they 
start to migrate immediately after introduction, as soon 
as they possess a sufficient amount of mass to exchange 
enough angular momentum with their surroundings to 
induce radial movement. Since the net torque acting on 
the companion is predominantly negative the direction 
of migration should be inwards in most cases. 

The companions are not allowed to accrete mass from 
the surrounding disk during the migration. However they 
may still trap disk material in their gravitational influ¬ 
ence radius (i.e. the companion’s circumplanetary disk), 
which might have an impact on the angular momentum 
transfer between the companion and the disk. Therefore 
while this additional mass is accounted for when comput¬ 
ing the torques responsible for gap formation, the addi¬ 
tional mass is not used to test the criterion directly. 

For each prepared SG disk set-up we introduce a 15, 
30, 45 and 80 Mjup companion separately. The simula¬ 
tions are performed four times for each companion rep¬ 
resented by different colors in Figure [2] Each time the 
companions are introduced with the same orbital radius 
but at different azimuthal angles (0, 7r/2, tt and 37r/2) to 
ensure the overall migration and gap-opening results are 
statistically significant. These simulations are run for 20 
orbits@ 

® This time period is a reasonable compromise between a short 
simulation time and acquiring sufficient data to distinctively char¬ 
acterize the companions’ migrations. 


In the MMSN-like disks we introduce companions with 
0.5 Mjup, 1 Mjup and 2 Mjup at 10 AU and allow them to 
migrate freely. The companions’ position is then tracked 
for 600 orbits (at 5 AU). 

4. RESULTS 

The most direct indication for gap-opening comes as a 
drastic slow-down in the migration rate. Also an order 
of magnitude decrease in the su rface density in the co¬ 
orbital region (|Crida et al.ll200^ can confirm gap clear¬ 
ing. Through this approach we assess whether or not 
the companions open gaps, and if so, at which orbital 
separation. 

4.1. SG disk simulations 

In the reference SG disk the migrations are very rapid. 
In all cases the companions reach the inner grid bound¬ 
ary in a few orbits. The 80 Mjup companion suffers a 
significant deceleration at an orbital separation of « 70 
AU. However the subsequent radial motion is still much 
more rapid than would be expected from a gap-opening 
procedure. Figure [2] shows the evolution of the orbital 
separation for the 30 Mjup and 80 Mjup companions in 
the two SG disk set-ups. Appendix |A] contains the or¬ 
bital evolution figures for the 15 Mjup and the 45 Mjup 
companions. 

In contrast to the reference SG disk, the lighter SG disk 
offers conditions which cause the companions to migrate 
much slower right from the point of introduction. They 
often undergo a further decelaration until they reach a 
quasi-constant orbit. No companions reach the inner grid 
boundary during the first 20 orbits of their migration 
(except the 15 Mjup companion once). Although they 
do not clear their horseshoe region to the extent of a full 
gap, the slow-down in the inward migration suggests that 
given enough time they would do so. This can be seen in 
Figure [H which depicts the surface mass density distri¬ 
bution around the embedded 30 Mjup companion in one 
of the lighter SG disk simulations. After 20 orbits the 
companion’s co-orbital region still appears to be in the 
process of being cleared out. To investigate this further 
we run this particular simulation longer (until 60 orbits), 
and find that after ^ 30 - 40 orbits the depth of the co¬ 
orbital region reaches a fully developed gap level. How¬ 
ever, the companion’s orbital separation still decreases 
significantly during the 60 orbits simulation suggesting 
that the companion might have not yet reached a type 
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Fig. 2. — Evolution of the orbital separation of a migrating 30 Mjup (top) and 80 Mjup (bottom) companion in the reference (left) and 
lighter (right) SG disks. The simulations are performed four times starting at the same orbital radius but each run at a different azimuth 
leading to four separate planetary trails (displayed with different colors). In the reference SG disk the companions migrate rapidly through 
the disk preventing any potential gap-opening. In the lighter SG disk they decelerate enough or even remain at their introductory radial 
location so that a gap might evolve. 


II migration regime, commonly associated with a fully 
developed gap. We discuss the possible reasons behind 
the slow migration rates in the lighter SG disk in section 

[Ql 

A snapshot of the 30 Mj^p companion in the reference 
and lighter SG disks after 3 orbits is provided by the 
upper panels of Fig. |6l The companion in the reference 
SG disk is located at a lower orbital radius due to its 
more rapid inward migration. 

We conclude that in the lighter SG disk our empirical 
study appears to be generally consiste nt with the torque 
balance criterion of iCrida et al.l (|2006D (see Fig. [H mid¬ 
dle panel). In the reference SG disk however, despite 
being expected to open gaps (Fig. (U left panel), the com¬ 
panions migrate inwards rapidly. To investigate whether 
the time aspect might be a limiting constraint we conduct 
fixed orbit simulations and determine the gap-opening 
timescale for our SG disk set-ups empirically, as analyti¬ 
cal estimates for th e gap-opening timescale are based on 
linear analysis (e.g. iLin fc PapaloizoullIQSfil their eq. 9) 
while our simulations are in the non-linear regime. We 
then compare these with the migration timescale for the 
’’free” companions. 



Fig. 3.— Evolution of the disk surface density perturbation due 
to a migrating 30 Mj^p companion in the lighter SG disk. A gap 
is believed to be established when the surface density drops by ap¬ 
proximately one order of magnitude in the companion’s horseshoe 
region, which is satisfied at approx. 30 - 40 orbits after introduc¬ 
tion but not yet after 20 orbits. During the run the companion’s 
orbital separation decreases only slowly due to a strong attenuation 
of torques in the cleared co-orbital region. 
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orbital radius at 70 AU (left) and 140 AU (right) in the reference SG disk. A gap is assumed to be fully established when the surface density 
drops by approximately one order of magnitude in the companion’s horseshoe region. For the embedded 30 Mjup (80 Mj^p) companion a 
gap is formed after approximately 8 (6) orbits at 70 AU and after 35 (20) orbits at 140 AU. The time it takes for a gap to form is mostly 
larger than the respective migration timescale shown in Figure [2](see also Table 0 for both companions. 


4.1.1. Companions on a fixed orbit 

Using the reference SG disk set-up we firstly employ as 
before 15 Mjup, 30 Mjup, 45 Mjup and 80 Mjup conpan- 
ions at fixed orbital separations of 70 AU and 140 AU. 
The results for the 30 Mjup and 80 Mjup companions 
are depicted in Fig. |4] (and for the curious reader the 
results for the 15 and 45 Mjup companions are presented 
in Appendix lAl). 

At both introduction locations all companions succeed 
eventually at opening gaps. This is consistent with the 
general criterion’s prediction (see Fig. [U left panel). The 
15 Mjup companion is special in the sense that it is not 
predicted to clear a gap at 140 AU but does so in our 
simulation. However, taking the whole perturbing mass 
~ 22.7 Mjup, i.e. the companion mass and the mass 
of its circumplane tary disk (materi al within « 2/3 i?H 
of the companion: iCrida et al.l[200^ . eliminates the dis¬ 
crepancy between simulation results and prediction as 
this value lies slightly above the mass limit shown in Fig. 
[TJ We hence find that when radial migration is disabled 
and timescales are not considered, gaps open as predicted 
by the torque balance criterion when otherwise they do 
not. This provides a first indication that the timescale 
for gap-opening may be important. 


Secondly, all of the companions require more time for 
gap-clearing than a ’’free” migration would take. We es¬ 
timate the gap-opening timescale to be the average time 
taken to clear the inner and onter horseshoe region such 
that the surface density drops by approximately one or¬ 
der of magnitude in t he gap region - con sistent with the 
definition of a gap in ICrida et all (1200611 . This is suffi¬ 
cient to provide a reasonable estimate of the gap-opening 
timescale. Fignre [4] shows that the 30 Afjup (80 Mjup) 
companion clears a gap after « 8 (6) orbits at 70 AU 
and after m 35 (20) orbits at 140 AU. The migration 
timescales from 140 AU to the inner grid at 20 AU for 
the 30 Mjup and 80 Mjup planets are ~ 5 - 6 orbits and ~ 
11-13 orbits, respectively, and so the crossing timescale 
across a gap region that a companion may try to form is 
much less than this. We define the crossing timescale to 
be 


\dR/dt\ ^ ’ 

where i?HS is the half-width of the companion’s horse¬ 
shoe region or equivalently the distance over which 
the companion would need to migrate to start the 
gap-opening process again. Although it is not clear 
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TABLE 4 

Comparison between the gap-opening timescale tgap versus the crossing timescale tcross empirically determined at 70 and 
140 AU. The time is normalized to orbits at 100 AU for the SG disks and to orbits at 5 AU for the MMSN-like disks. Since 
IN the vast majority of the simulations tgap 3> tcross CAP-OPENING MIGHT BE PREV ENTED EVEN T H OUGH THE COMPANIONS ARE 
PREDICTED TO OPEN GAPS BY THE TORQUE BALANCE CRITERION OF ICRIDA ET AL.I II2006I1 . 


disk model 

planet mass [Mj^p] 

orbital separation [AU] 

fgap [orbits] 

tcross [orbits] 

reference SG 

15 Mjup 

70 

~ 12 

~ 1 

140 

~ 65 

~ 2 

30 Mjup 

70 

rw 8 

~ 1.2 

140 

~ 35 

~ 3 

45 Affjup 

70 

rw 8 

~ 1.5 

140 

~ 20 

~ 2 

80 Mjup 

70 

6 

8 

140 

~ 20 

~ 5 

heavier MMSN-like 

0.5 Mjup 

5 

400 

~ 3 

10 

> 600 

~ 30 

1 Mjup 

5 

~ 200 

~ 2 

10 

> 600 

~ 20 

2 Mjup 

5 

60 

~ 3 

10 

400 

~ 15 


TABLE 5 

Comparison between the empirically determined gap-opening timescale tgap for the 30 Mjup companion in the reference 
AND LIGHTER SG DISKS (tOp) AND FOR THE 2 Mjup COMPANION IN THE HEAVIER AND LIGHTER MMSN-LIKE DISKS (bOTTOm). ThE TIME IS 
NORMALIZED TO ORBITS AT 100 AU FOR THE SG DISKS AND TO ORBITS AT 5 AU FOR THE MMSN-LIKE DISKS. ThE GAP-OPENING 

TIMESCALE APPEARS TO BE INDEPENDENT OF DISK MASS. 


planet mass [Mjup] 

orbital separation [AU] 

tgap [orbits] 

reference SG disk 

tgap [orbits] 
lighter SG disk 

30 Mjup 

70 

8 

~ 8 

140 

~ 35 

~ 35 

planet mass [Mjup] 

orbital separation [AU] 

tgap [orbits] 

heavier MMSN-like disk 

tgap [orbits] 

lighter MMSN-like disk 

2 Mjup 

5 

60 

~ 60 

10 

400 

~ 450 


what a reasonable estimate for flns is in our mass 
regi me, we extrapolate the results o fIMasset et alJ (|2006f) 
and iPaardekooner fc Papaloizonl (j2009fl . which found 
through streamline analysis that i?HS ~ 2.5 i?H for higher 
mass planetfl The migration rate dR/dt is defined as 
the gradient of the migration curves in Figures [H [5] and 
lA.fl at the relevant location. Table IH compares the gap¬ 
opening timescale (for planets on a fixed orbit) to the 
crossing timescale determined from the migration simu¬ 
lations at particular disk radii and shows that the sim¬ 
ulations in which gaps do not open are those in which 
the crossing is too fast. Since in most cases tcross < tgap, 
this may explain the absence of gap-opening in the free 
migration simulations. 

^ IMasset et al.l II2006I) fi tted data points up t o a planetary mass of 
0.2 Mt„p. IPaard^ooper &: Papaloizonl II2009I1 used the dimension¬ 
less parameter q/h^ to represent their planet masses and examined 
a mass regime below q/h^ = 5. As a comparison: our MMSN-like 
simulati ons correspond to values q/h^ = 1.4 — 5.7, which is almost 
fully in IPaardekooper fc Papaloizonl Il2009ll ’s examined regime. 
The SG disk simulations are in the range q/h^ = 6.8 — 224 and 
hence involve a large extrapolation of IPaardekooper &: Papaloizonl 
Il2009l )’s results. 


As an aside we also point out that the gap-opening 
timescale decreases with increasing companion mass and 
decreasing orbital separation. This is consistent with 
the prediction of the torque balance criterion, which is 
satisfied more easily by massive companions at smaller 
orbital radii (see Fig. [T]). Furthermore we observe no 
significant deviation in the gap-opening timescale when 
changing the disk mass even in the non-linear regime 
- consistent with the fin dings of the linear analysis by 
iLin fc Papaloizoul (Il986f) (s ee T able [5] and compare the 
top panels of Figures HI and IA.41) . 

4.2. MMSN-like disk simulations 

Figure [5] shows the evolution of the companions’ or¬ 
bital separation in the lighter (left), intermediate (mid¬ 
dle) and heavier MMSN-like disks (right). As shown in 
Fig. [T] (right panel) the torque balance criterion predicts 
that the 2 Mj^p companion should immediately open a 
gap, the 1 Mjup companion should do so somewhat fur¬ 
ther in and the 0.5 Mj^p companion should do so after 
having crossed large parts of the disk. In the lighter and 
intermediate disks the 1 Mjup and 2 Mjup companions 
confirm this prediction. They both start to slow their 
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Lighter MMSN-like disk 
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Intermediate MMSN-like disk 



Heavier MMSN-like disk 



Fig. 5.— Evolution of the orbital separation of a migrating 0.5 Mjup (solid line), 1 Mjup (dotted line) and 2 Mjup companion (dashed 
line) in the lighter (left), intermediate (middle) and heavier MMSN-like disks (right). The larger the disk mass the more rapid the inward 
migration of the companions. In the lighter and intermediate disks we observe the predicted gap-opening by the 2 Mj^p companion whereas 
the 1 Mjup companion only partially manages to clear its co-orbital region during the simulated 600 orbits. The 0.5 Mjup companion shows 
little signs of gap clearing due to its rapid inward migr ation in all set-ups. The companions’ ability to open a gap depends significantly on 
the disk mass, although the torque balance criterion of ICrida et al.l 11200611 does not suggest this. 



-10 -5 0 5 10 -10 -5 0 5 10 

X axis [AU] X axis [AU] 

Fig. 6. — A snapshot of the inwards migration of a companion through the disk. Top left: 30 Afj^p companion in the reference SG disk 
after 3 orbits (at 100 AU). Top right: 30 Mjup companion in the lighter SG disk after 3 orbits (at 100 AU). Bottom left: 2 Mjup 
companion in the lighter MMSN-like disk after 15 orbits (at 5 AU). Bottom right: 2 Mjup companion in the heavier MMSN disk after 
15 orbits (at 5 AU). The lack of disk material in the co-orbital region is visible by the darker colors around the embedded companions. 
The companions in the reference SG and heavier MMSN-like disks are located further in than their counterparts in the lighter SG and 
MMSN-like disks since they migrate inwards more rapidly due to a stronger torque (see Figure [7|l. 
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migration by clearing their co-orbital region early. How¬ 
ever, the 0.5 Mjup migrates to the inner grid boundary 
and indicates no gap clearing. In the heavier MMSN-like 
disk, the migration times of all the companions decrease 
substantially. The bottom two panels of Figure [S] depict 
the 2 Mjup companion at 15 orbits after introduction 
into the lighter (left) and the heavier disks (right). The 
2 Mjup companion is located at a smaller orbital radius 
in the heavier MMSN-like disk compared to the lighter 
MMSN-like disk due to its more rapid migration. 

As with the high-mass SG simulations (see section IQ) 
we find that the disk mass is a non-negligible factor for 
migration and also affects the companion’s ability to 
open a gap. This is surprising as the disk mass is not 
includ ed in the torque balance criterion of iCrida et al.l 
(j200(itl and thus its prediction is the same for the differ¬ 
ent MMSN-like disks. 

4.2.1. Companions on a fixed orbit 

As with the SG disks in Section 14. 1.11 we perform sim¬ 
ulations with the companions held on a fixed orbit at 10 
AU and 5 AU and find that the gap-opening timescale 
is much larger than the migration timescale in the heav¬ 
ier MMSN-like disk (see Tabled]). This confirms that the 
timescale argument is the limiting factor for gap-opening. 
We also determine the gap-opening timescale for the 2 
Afjup companion in the lighter MMSN-like disk and find 
that the timescale is independent of the disk mass (see 
Table O as with the SG disks simulations. The curi¬ 
ous reader finds the according surface density evolutions 
shown in Figures [A.31 and lA.41 

4.3. Migration Analysis 

We have found that companions in disks massive 
enough migrate very rapidly inwards even if they are 
predicted to open gaps and slow down thei r migration, 
accor ding to the torque balance criterion (|Grida et al.l 
1200611 . The results in Table m show the reason behind 
this discrepancy: the companions simply lack the time 
to open gaps, as the gap-opening timescale is larger than 
the crossing timescale (eq. fT^ . We find that in both the 
SG and the MMSN-like disks the torque balance con¬ 
dition for gap-opening does not seem to sufficiently de¬ 
scribe when gap-opening will occur. How can such rapid 
migrations in some of our disk models be explained? 

4.3.1. Torque 

To provide an explanation for the very rapid migra¬ 
tion timescales we investigate the tidal torque between 
the companion and disk material. We evaluate the 
torque from each grid ce ll on the companion given by 
(|de Val-Borro et at 1200611 

T = RcGMcMceii-^— ^ , sm{Rc,d). (14) 

+ Cc) ' 

Here d represents the distance between the companion 
and the cell, sm{Rc,d) is the sine of the angle between 
the Rc and d vectors and Cc = 0.7H (as for the computa¬ 
tion of the companion’s gravitational potential). FigurejTj 
displays the azimuthally summed up contributions of the 
torque acting on the companion at a particular orbital 
radius in our migration simulations of the 30 Mj^p com¬ 
panion in the reference (top left) and lighter (top right) 


SG disks after 3 orbits, as well as the 2 Mj^p companion 
in the lighter (bottom left) and heavier (bottom right) 
MMSN-like disks after 15 orbits. To examine the impact 
of the circumplanetary disk on the torque we also show 
the torque profile with this region’s exclusion in red. We 
find in general that the corotation torque, exerted by disk 
material within the horseshoe region, plays an important 
role on the magnitude and the orientation of the total 
torque. In the reference SG and heavier MMSN-like disks 
the torque is predominantly negative and the circum¬ 
planetary disk has also an additionally negative effect 
(compare red and black profiles in Fig. [7|). Such a strong 
negative torque would explain the rapid inwards migra¬ 
tions and very short migration times in these simulation 
runs. Torque profiles with a large peak in the immediate 
proximity of the com panion may be indicative of a Type 
III migration regime (jMasset fc PaDaloizoiJ200^ . In the 
lighter SG disk the torque is generally weaker than in the 
reference SG disk. The circumplanetary disk has also a 
strong effect but in this case provides both positive and 
negative contributions, which may cancel each other out 
and thus provide a possible explanation for the slower mi¬ 
gration of the 30 Mjup companion in the lighter SG disk. 
In the lighter MMSN-like disk the torque is substantially 
weaker than in the heavier MMSN-like disk. Moreover, 
the torque is more equally distributed and hence the ex¬ 
clusion of the circumplanetary region has a signihcantly 
smaller effect on the torque profile, which is consistent 
with a common Type I migration regime. 

Although not depicted here, examining the tidal torque 
profile of the 0.5 Mjup companion at several stages of its 
inward migration in all three MMSN-like disks suggests 
a Type I to Type HI regime change during the inward 
migration. This is visible in Fig. [5] by the companion’s 
acceleration of the migration rate between an orbital sep¬ 
aration of 6 - 9 AU (depending on the disk model). 

4.3.2. Mass deficit and timescales 

There are two dynamical effects that play a dominant 
role in the corotation torque of a migrating companion: 
A rapid inward migration of a companion causes disk 
material, which otherwise would be trapped on a horse¬ 
shoe orbit, to flow across the companion’s corotation re¬ 
gion. The associated negative angular momentum trans¬ 
fer from the crossing particles onto the companion scales 
with the drift rate. In contrast, material within the com¬ 
panion’s circumplanetary region, dragged inwards with 
the companion, yields a positive torque on the compan¬ 
ion. Hence this mechanism slows down the migration. 
From these considerations it is plausible, that the mass 
found within the horseshoe region as well as the mass 
trapped in th e circumplanetary regi o n imp acts the mi¬ 
gration rate. iMasset fc Panaloizoij (j2003ll stated that 
for a very rapid migration (Type HI migration) to estab¬ 
lish ,the ” co-orbital mass deficit ’O, Sm, should be of the 
same order of magnitude as the companion’s mass plus 
the mass within the circumplanetary disk. Me + Mcp. 
Moreover, if Sm > Mc + Mcp the positive feedback of the 
crossing disk material causes the migration to become ex¬ 
tremely fast, a regime known as ’’runaway” migration. If 
Sm < Mc + Mcp the negative feedback of the circumplan- 

The reduction in disk mass in the companion’s horseshoe re¬ 
gion (with width 5 Rh) due to its presence 
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Orbital separation [AU] Orbital separation [AU] 

Fig. 7. — Radial distribution of the total tidal torque acting on the companion due to the influence of the disk material. The effect of 
excluding the torque from the circumplanetary disk is shown in red. Vertical lines depict the companion’s position (dot-dashed line) and the 
theoretical outer and inner horseshoe region boundaries at di 2.5 Hill radii (dotted lines). Top left: 30 Afj^p companion in the reference 
SG disk after 3 orbits (at 100 AU). Top right: 30 Mjup companion in the lighter SG disk after 3 orbits (at 100 AU). Bottom left: 2 
Afjup companion in the lighter MMSN-like disk after 15 orbits (at 5 AU). Bottom right: 2 Mj^p companion in the heavier MMSN-like 
disk after 15 orbits (at 5 AU). In the reference SG and heavier MMSN-like disks the corotation torque is predominantly negative and 
the circumplanetary disk also provides a negative torque contribution. In the lighter SG disk the torque is generally weaker and the 
circumplanetary disk has a strong effect on the torque profile, providing both a positive and a negative torque. In the lighter MMSN-like 
disk the torque is weak, more equally distributed and the circumplanetary disk has only a small impact on the torque profile. This profile 
is consistent with a common Type I migration regime. These four graphs correspond to the same images displayed in Figure [6] 


TABLE 6 

Companion mass including the mass of its circumplanetary disk, Me + Mep, co-or bital mass deficit, Sm, libration timescale, 

riib, AND MIGRATION TIMESCALE, Tmigr, FOR THE SIMULATIONS DESCRIBED IN SECT. 14.3.11 We FIND THAT Sm IS SOMEWHAT LOWER THAN 
Me + Mep IN ALL THE SIMULATION RUNS BUT STILL OF THE SAME ORDER OF MAGNITUD E. ThI S INDICATES THAT THE MIGRATIONS MIGHT BE 
WITHIN THE GENERAL TyPE III REGIME BUT NOT IN THE RUNAWAY REGIME (SEE SECT. 14.3.21) . In CONTRAST, THE MIGRATION TIME Tmigr 
CORRESPONDS FAIRLY WELL TO THE LIBRATION TIMESCALE Tub IN THE REFERENCE SG AND HEAVIER MMSN-LIKE DISKS SUGGESTING INDEED 

A Type III runaway migration. 


simulation run 

Me Mep 

5m 

T-migr [orbits] 

Tiib [orbits] 

30 Mjup, reference SG disk, 3 orbits 

~ 43.8Mjup 

~ 17Mjup 

~ 5 

4 

30 Mjup, lighter SG disk, 3 orbits 

~ 37.5Mjup 

~ 17Mjup 

> 60 

4 

2 Mjup, heavier MMSN-like disk, 15 orbits 

~ 3.2Mjup 

~ 2.1Mjup 

~ 24 

~ 17 

2 Mjup, lighter MMSN-like disk, 15 orbits 

~ 2.07Mjup 

~ 0.7Mj,p 

> 600 

~ 17 
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etary disk prevails and weakens the negative corotation 
torque enough to prevent any runaway. The migration 
timescale in the runaway case is thought to be a few times 
the libration timescale of dis k particles on a horsesh oe 
orbit around the companion (jLin fc PaDaloizoull20ir)tl . 

We perform a rough analysis of the mass deficit and 
the migration versus libration timescale for the simula¬ 
tions in question to determine whether our simulations 
fit the above stated conditions for a Type III migration 
regime in general or even specifically a runaway migra¬ 
tion. First, we compute the mass deficit by 


Rc+Rhs 


Sm = 27 t 


AERdR, 


Rc—Rns 


(15) 


where AS is the difference between the unperturbed 
surface density (i.e. the ’’smoothed” density with the 
stochastic peaks and troughs flattened out at the time 
of the companion’s introduction) and the actual surface 
density at a given rB. Second, we wr ite the libration 
timescale as (iBaruteau fc Mas^l2008all 


SttRp 

3npi?HS 


(16) 


In addition, we estimate M^p, the mass of the circum- 
planetary disk. 

We focu s our attention on the simulations described 
in section 14.3.11 and shown in Fig. [6] and Fig. 0 For 
these we compare the companion’s mass plus the mass 
of the circumplanetary disk, the mass deficit, the libra¬ 
tion timescales, and the migration timescales in Table [HI 
The values for Me + M^p and 6m are the average be¬ 
tween the 2nd and the 4th orbit in the massive SG disk 
simulations and between the 11th and the 20th orbit in 
the MMSN-like disk simulations. The mass deficit is also 
visible in Figure [6] (darker regions around the embedded 
companions). 

The general condition for Type III, i.e. Sm ~ Me -I- 
Mep, is satisfied roughly for all the depicted cases. Al¬ 
though not strongly, it is noticeable that the difference 
between Sm and Me + Mep is smaller in the reference SG 
disk and heavier MMSN-like disk compared to the lighter 
SG disk and lighter MMSN-like disk, respectively. The 
condition for runaway migration, i.e. Sm > Me + Mep, 
is not satisfied for any of our simulations. In contrast to 
the condition imposed on the mass deficit, the timescale 
comparison suggests that the reference SG and heavier 
MMSN-like disks might be in the runaway regime as their 
migration times are comparable to the respective libra¬ 
tion timescale. Overall it remains unclear if the very 
rapid migration regime observed in many of our runs be- 


Our method of calculating the mass deficit is similar, 
but slightly differen t to that often used in the literature (e.g. 
IKlev &: Nelsonl2012IV The latter assumes a uniform surface density 
over the horseshoe region equivalent to that at the inner separa- 
trix. We find that while this approach is accurate for light disks 
and small planets, it largely overestimates the total mass the horse¬ 
shoe region would have without the influence of the companion in 
our simulation set-up. For example, for the 30 Mjup companion in 
our SG disks, this equates to the mass deficit being almost equiv¬ 
alent to the entire disk mass. Therefore a more accurate way of 
calculating the mass deficit in this region is needed rather than 
assuming a uniform surface density. 


Heovier MMSN-like disk 



Time [Orbits at 5 AU] 

Fig. 8.— Evolution of the orbital separation of a migrating 0.5 
Mjup (solid line), 1 Mj^p (dotted line) and 2 Mj^p companion 
(dashed line) in the heavier MMSN-like disk with different soft¬ 
ening lengths. Black: ec = 0.7, esq = 0.8. Red: Ec = 0.3, 
^SG = 0.4. A lower softening length causes a more rapid migra¬ 
tion. 


longs to the categories referred to as ” runaway” or simply 
’’Type III migration” in the literature. 

4.4. Numerical tests 
4.4.1. Resolution 

Since rapid migration depends on the local torques 
near the corotation region, inadequate resolution can 
affect the migration results. We repeat the simulation 
involving the 2 Mjup mass companion in the heavier 
MMSN-like disk, as well as the 30 Mjup companion in 
the reference SG disk using twice as many grid cells in 
both the radial and azimuthal directions and find that 
this does not affect the migration timescales. 

4.4.2. Softening length 

In their Fig. 7, iMiiller et al.l (j2012t) investigated the 
impact of Cc on the net torque acting on their planet. 
Varying Cc from 0.4i7 to l.OM, they find that the torque 
density in the planet’s influence region decreases by a 
factor ^ 2 to 3 for a Neptune mass planet in a disk with 
scale height, H = 0.05i?. 

Since the net torque on the companion is the ’’motor” 
of migration we perform softening length tests by con¬ 
ducting another suite of simulations using the heavier 
MMSN disk set-up as in section but with Cc = 0.3 
and esG = 0.4. We find that a lower softening length 
increases the angular momentum exchange between the 
companion and disk by diminishing the smoothening out 
in the companion’s horseshoe region. This causes the 
companions to migrate even more rapidly through the 
disk (see Fig. [5]). 

Our results are b y such consistent with the findings of 
IMiiller et al.l (|20I2D . This shows the importance of using 
an accur ate softening l ength in migration simulations. 
Based on IMiiller et al.l (j2012t) we choose Cc = 0.777 and 
esG = 0.877 since it provides a reasonable description of 
the gravity force. 
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4.4.3. Growth time after companion introduction 

We introduce the companions by linearly increasing 
their mass over one orbit (see section ESI). To establish 
the effect of this numerical choice we repeat the simula¬ 
tions with the 30 Mjup companion in the reference SG 
disk and the 2 Mj^p companion in the heavier MMSN- 
like disk with introductory growth times of 10 and 600 
orbits, respectively. We find that as soon as the com¬ 
panion accumulates enough mass to establish an effective 
torque on the surrounding disk the subsequent migration 
is as rapid as in our standard runs. Hence we conclude 
that the salient features of our simulation results are in¬ 
dependent of the numerical growth time. 

5. DISCUSSION 

The gap-opening criterion is fundamental to under¬ 
stand planet survival from theoretical and observational 
perspectives. We find that not only is it important to 
consider whether a gap can open in a disk, but more 
importantly whether it can open qui ckly enough. The 
frequently used gap-opening criterion (jCrida et al.ll2006[l 
was derived for a planet on a fixed orbit and not for 
a migrating planet. Our results show that even when 
the gap-opening criterion is satisfied there are still in¬ 
stances when a gap does not open. The timescale on 
which the planet carves out the gap, tgap, is the gap¬ 
opening timescale. However, at the same time, the planet 
is also migrating through the region that it is carving 
out on a timescale, Uross- If the planet crosses this 
region quicker than the time it takes to carve out the 
gap, i.e. if tcross < tgap, a gap will not form. Indeed, 
when comparing the gap-opening timescales to the mi¬ 
gration timescales in our simulations, we find that in 
many cases even when th e gap-opening criterion accord¬ 
ing to ICrida et'al] (|2006[1 is satisfied the gap simply can¬ 
not open because the timescale constraint is not satisfied. 
Our results show that both a torque balance criterion and 
a timescale condition are necessary. The combination of 
the two will always catch the limiting factor, resulting in 
an equivalent or stricter criterion for gap-opening. Thus 
previous assessments of the gap-openi ng criterion where 
the planet is kept on a fixed orbit (e.g. ICrida et'ani2006L 
iDuffell &: MacFadvenll2013ll may not be sufficient enough 
to determine if a gap will form in a disk. We also find 
that gap-opening is dependent on the disk mass when 
this is not a factor in the semi-analytical torque balance 
equation (but is implicit in the timescale condition). 

We can understand from a phenomenological perspec¬ 
tive why this additional timescale condition is neces¬ 
sary to consider. The crossing timescale depends on the 
torques acting on the planet which depend on the planet- 
to-primary mass ratio, the disk location, aspect ratio, the 
surface mass density and hence disk mass, the surface 
mass density profile and the temperature profile. The 
gap-opening timescale depends on the aspect ratio, the 
planet-to -primary mass ratio, the viscosity and th e radial 
loc^ i on llHouriga n fc Wa rd 11984 IRafikovI l2002t lEdead 
I2007t iBaruteau et al.l 120131). Out of these py ameters, 
the current gap criterion bv ICrida et al.l (|2006ll does not 
consider the disk mass or the surface mass density and 
temperature profiles, though intuitively they would be 
expected to affect gap-opening. Yet these are crucial to 
determining a planet’s migration rate through the region 


that it is trying to open a gap in. In fact the migration 
and gap-ope ning timesca l es are affected by all the vari¬ 
ables in the ICrida et al.l (|2006f) criterion as well as the 
aforementioned additional variables. 

Furthermore, our results have a broad impact since we 
show that this timescale condition is important to con¬ 
sider in both high mass self-gravitating disks and low 
mass T Tauri disks where planets/companions are ex¬ 
pected to form by gravitational instability and core ac¬ 
cretion, respectively. 

Note that a timescale criterion for gap-opening of the 
form tcross > tgap would need to naturally take into ac¬ 
count very rapid migration (e.g. Type HI). This occurs 
for intermediate mass planets or for planets in very tur¬ 
bulent disks which end up having circumplanetary mate¬ 
rial around them that causes the corotation torque to be 
more important. Different aspects of our analysis suggest 
that the planet evolves in a regime which exhibits some 
of the features attribute d to Type Hl-like migration , or 
even runaway migration (iMasset fc PaDaloizoull2003ll . In 
these cases the corotation torques yield the most impor¬ 
tant contribution and govern planet migration. The rel¬ 
atively large effective viscosity arising in gravitationally 
unstable disks may explain why the corotation torques 
never saturate in this case, leading to Type HI migra¬ 
tion: mass may simply be transferred to the co-orbital 
region faster than the libration timescale. This in turn 
prevents gap-opening, not because the gap is refilled by 
the viscous flow of material through it but because the 
planet transits too rapidly through the region where it is 
trying to open a gap. 

We note that disks that have different structures will 
affect a planet’s migration rate and direction, and thus its 
ability to open a gap. Furthermore, interactions between 
multiple companions, as well as the interactions between 
companions and other disk structures such as the spi¬ 
ral structures caused by the presence of further compan- 
ion(s) or by self-gravity in high mass disks will affect the 
torques, possibly causing fluctuations on short timescales 
and adding a significant non-linearity to the outcome. 
These may well cause gap-opening to be harder. Even 
Type I migration might not display a smooth continu¬ 
ous radial motion throughout the disk, but may consist 
of periods with increased or decreased migration rates 
due to inhomogeneities in the disk structure. A cap- 
tu re of a planet into a z ero torque location as described 
bv iBitsch fc KlevI (|2011h or into mean motion resonance 
with another planet may well aid the gap-opening process 
due to the long migration timescales in such scenarios. 
In all of these cases, the timescale criterion is still a rele¬ 
vant condition to satisfy since these can be reassessed in 
the context of differing situations occurring in the same 
background disk. 

5.1. Comparisons with previous work 

The idea of a timescale constraint has been considered 
previously. In this paper we stress its applicability to 
a wide range of disks including young massive gravita¬ 
tionally unstable disks as well as older low mass T Tauri 
disks, i.e. disks in which planets may form by gravita- 
ti onal instability or c or e accr etion. 

iHourigan fc Wai^ l)1984ll and iWard fc HouriganI 
(|I989ll suggested that for a gap of a particular width 
to open, the planet must drift across that region more 
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slowly than th e time it takes t o open the gap, i.e. 
Across > ^gap- iLin fc Papaloizoul (|1986ll on the other 
hand suggested that since the presence of a planet 
causes the surface density in the disk to vary the 
relevant timescales to consider are the crossing timescale 
and the timescale to modify the surface density, given 
by = (H/R)tgap, and that a gap would form if 
Across > tA - This is an easier criter ion for gap-opening 
than what iHourigan &: Ward! (|1984ll suggested. Given 
that the migration simulations performed in this work 
are in the non-linear regime, the analytical expressions 
provided by these authors cannot be used to directly 
compare with the empirically-obtained timescales in 
this paper, nor do our results indicate exactly which 
timescale condition should be used. We simply point out 
the importance of considering a timescale criterion in 
additi on to the torque balance criterion of ICrida et al.l 

(I2OOI 

iTakeuchi et al.l (I1996I1 also explore the gap-opening cri¬ 
terion and show that gaps are formed when two condi¬ 
tions are satisfied (see their Figure 8): (i) tcross > igap 
and (ii) topen < ^ciose- The latter criterion is essentially 
the same as the criterion that the viscous torques must 
balance the gravitational torques. They also show that 
when the latter criterion is satisfied, gaps do not neces¬ 
sarily open unless the former is satisfied. Our results are 
therefore in ag reement with theirs . 

Sinr ilar to IHourigan fc Ward! (|1984ll . IMuto et'H] 
(|2010ll also discuss the idea that a gap must open only 
if the migration timescale is longer than the gap-opening 
timescale. Using the Type I migration timescale and 
their numerical results for the gap-opening timescale us¬ 
ing local 2D shearing sheet simulations, they find that 
it is easier for planets to open gaps, i.e. lower mass 
planets are able t o open gaps compared to what the 
ICrida et'all (|2006ll criterion predicts, in contrast to our 
results. However, it must be noted that these were for 
inviscid disks. Since the viscosity will affect both the 
gap-opening timescale (lengthening it) and the migra¬ 
tion times cale (higher v iscosities decrease the migration 
timescale; lEdgad 1200711 . the inclusion of viscosity will 
make the timescal e criterion harder to satisfy than what 
IMuto et aP (|2010ll suggest for inviscid disks (though we 
note that they do specify that the inclusion of viscos¬ 
ity does require an additional condition - also based on 
timescale arguments - to be satisfie d). Note that the 
Type I migration timescale used by IMuto etall 1)201011 
was an appropriate choice for them since they consider 
low mass (approx. 1-23 M 0 ) planets but is not an ap¬ 
propriate choice for our simulations since the migration 
o f our high mass co mpa nions is in the no n-linear regime. 

IZhu et aLKlMTl and IVorobvovI (|2013ll perform hydro- 
dynamical simulations of fragments in gravitationally un¬ 
stable disks and find that in some cases they can open a 
gap and survive. Both authors showed that the torque 
balance criterion is satisfiecQ- It would be particularly 
interesting to determine whether the timescale condition 
in their simulations is also satisfied. 

Note that t he authors use different torque balance criteria: 
|Zhu_e^_aD 120121) requir e that Rff > H and Me > AOMtOib? 
IILin fc P aDaloizou l 119931 1 while IVorobvovI 1120131 use the criterion 
bv lCrid^^^am2006l) except that they use the mass of the fragment 
plus the mass within a Hill radius (i.e. that of the circumplanetary 
disk). 


iZhang et al.l l)2014[l examined how self-gravity and disk 
mass influence the ability of a Jovian-mass planet to 
open gaps. They used disks in the mass range of the 
MMSN model. They found that disks that are too mas¬ 
sive may inhibit gap-opening, consistent with our re¬ 
sults. Whereas we attribute this to a timescale problem 
they argued using ID hydrodynamical simulations that 
the gravitoturbulence in a disk with a surface density 
E > 3.5MMSN suppresses gap formation. In 2D simu¬ 
lations their planet migration time scales inversely to the 
disk density. In this sense our results are consistent with 
this trend. 

5.2. Additional considerations 

We note that in the results presented here the planet 
has a fixed mass. Migration, in any regime, and gap 
opening both depend on planet mass. As a planet grows, 
it may pass through a rapid migration phase before slow¬ 
ing down as it starts to open a gap. Accretion onto 
the planet (and hence removal of gas in the surround¬ 
ing r egion) can promote gap formation (jBrvden et al.l 
II999H which will also have an effect on the gap-opening 
timescale and hence on whether a planet can open a 
gap quickly enough. On the other hand, although the 
planet’s mass affects the migration timescale, if the 
growth timescale is too long compared to the migration 
timescale, a planet’s growth may no t be affected (analo - 
gous to the black hole binary regime: ITvanov et '^11199911 . 
Regardless, the timescale criterion presented here would 
still be relevant in the case of planet growth, but it sim¬ 
ply means the the criterion woul d have to be re-a s sessed 
during the evolution. Recently, iMorbidelli et al.l (|2014fl 
confirmed earlier studies by finding a substantial mass 
accretion mechanism onto a planet on a fixed orbit sur¬ 
rounded by a gap. With a 3D model they found the gas 
flow in the gap varied significantly with vertical height 
above the mid-plane. Although their result does not de¬ 
scribe accretion onto a migrating planet, this suggests 
that in future studies the use of a 3D model might be 
ultimately necessary to realistically simulate the torque 
in the planet’s horseshoe region and hence gap-opening. 

We also note that in the SG disks the turbulence al- 
ph a parame t er, wh ile close to a = 0.013 as predicted 
by iGammid (|2001fl . may not remain at that value in 
the whole disk. If a varies then the prediction for gap¬ 
opening also changes. Figure [9] shows how the prediction 
changes in the reference self-gravitating disk if the turbu- 
len ce varies by a fa ctor of 2 compared to the expectation 
bv IGammid ((2001!) . While this does change the predic¬ 
tion, the companions should still open gaps during their 
inwards migration. Thus our conclusion that the gaps 
have not opened because of the timescale argument is 
still valid. 

5.3. Implications for planet formation models 

Our results may have important implications on the 
use of the criterion in studies using population syn¬ 
thesis models for planets formed by core accretion or 
gravitational instability, as well as the tidal downsiz- 
ing model for planets formed by gr avitational instability 
(|Bolev et al.ll2010l : lNavakshinll201C)!) . Gonventional popu¬ 
lation synthesis models for standard core accretion plan¬ 
ets typically assume that a planet migrates on a Type 
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Gop-opening prediction, reference SG disk 



Fig. 9. — Companion’s m ass versus or b ital ra dius indicating the 
torque balance criterion of ICrida et al.l II2Q06I) in red for the ref¬ 
erence SG disk using turbulence q parameters of 0.0065 (dotted), 
0.013 (solid) and 0.026 (dashed). Companions are expected to open 
a gap when situated above the limit. The crosses indicate where 
in the disk the companions are introduced in our simulations while 
the dashed lines indicate the inward migration tracks. The crite¬ 
rion is stricter if the turbulence parameter is increased but in all 
the simulations the companions are expected to open gaps at some 
stage during their inward migration. 


I migration timescale until it has grown large enough 
to satisfy the gap-opening criterion. At this stage the 
models then instantaneously make the pl anets migrate 
on a Type II migrat ion timescale (e.g. lAlibert et al.l 
[2(MlIda fc Lh][2f)0l . or more recently involve a smooth 
transition functi on between the Type I and Type II mi¬ 
gration regimes (|Dittkrist et al.l[2014ll . This is also the 
case for N-body simula tions of planetary systern forma¬ 
tion in evolving disks (|Coleman &: NelsCTl20I4| l. Pop¬ 
ulation synthesis for disk instabi hty planets has only 
recently begun to be d eyeloped (jForgan fc Ricd I20I31 
IGalyagni fc May^l20I4D . Both authors instantly change 
from Type I to Type II migration when a gap is deep 
enough to haye opened. The former assume this to be 
the case when the pressure stability condition (eq. [T|) 
is satisfie d^ while the latter assume that gap-opening 
occurs when iCrida et al.l l)2006ll ’s criterion (eq. [3|) is sat¬ 
isfied. 

Note that population synthesis models do not consider 
yery rapid migration scenarios like Type III, though in¬ 
corporating a gap-opening timescale condition into such 
models could naturally account for this, as it would natu¬ 
rally ensure that gap-opening would take longer than the 
typical migration timescale, thus providing a more accu¬ 
rate estimate of the survivability of planets. Our results 
show that gap-opening is the same or harder than that 
suggested by the torque balance criterion and therefore 
such population synthesis models are likely to underesti¬ 
mate how quickly a planet transitions from the Type I to 
Type II migration regimes, and thus have an important 
effect on a planet’s survivability. 

From the gra v itation al instability perspective, 
lOalvagni fc MavCT (j20I4ll found that the interplay 
between the collapse timescale, the migration timescale 

Note that their pressure stability condition marginally differs 
from equation Q in that the right hand side of their equation is 
I?p(5/2)(i/3). 


in the Type I regime, the gap-opening timescale and the 
tidal mass l oss timescale based on the tidal do wnsizing 
hypothesis (iBolev et al.ll20I0l : lNavakshinll20I0ll is what 
determines the fraction of survi ving gas giants . The y 
used the simulation results by iGalvagni et al.l (j20I2ll . 
who showed that the collapse to planetary densities 
occurs on timescales of ^ 10'* yr (comparable to Type 
I migration timescal es), and simply ad opted the torque 
balance criterion of iCrida et al.l (I2006D to decide when 
planets open a gap and enter the slow Type II migration 
regime. By comparing models with and without gap¬ 
opening they found very large qualitative differences in 
the fraction of surviving planets, mass distributions and 
semi-major axes after less than a million years of disk 
evolution. In particular, without gap-opening < 10% of 
the initial population of gas giants survived, and most 
of them ended up on close-in orbits consistent with 
Hot Jupiters, whereas with gap-opening the fraction of 
surviving planets increased to > 50% , with a sizeable 
fraction on orbits having R > 5 AU. The final mass 
distribution also differed, being skewed to lower masses 
when gap-opening is not included as tidal mass loss 
is more effective when planets migrate inward fast. 
Other semi-analytical models studying the evolution of a 
population of disk instability planets, incl uding also dust 
and p lanetesimals, were published by iForgan fc Ricd 
(|20I3l) . They considered an even simpler mass threshold 
above which gap opening occurs depending on the disk 
pressure scale height and adopted standard Type I and 
Type II migration timescales based on the results of Bate 
et al. (2003). Our results show that gap-opening should 
occur on timescales shorter than 10® yr, corresponding 
to the typical clump collapse timescale, in order to 
affect a planet’s orbital evolution, and that this requires 
relatively low mass disks (Mdisk < 0.05Mq) for planets 
with masses of a fe w Mj nn. M odels such as t h ose o f 
IGalvagni fc Mave'^ (j20I4l ) and iForgan fc Bicel (j20I.3f l 
should thus be revisited by introducing a criterion for 
the gap-opening timescale, and possibly mechanisms of 
disk dispersal acting since we show that the disk mass 
clearly plays a major role. 

5.4. Implications for observations 

There are many observations that indicate the presence 
of substantial gaps in gas-rich disks of the order of 10 
AU or larger. These gaps have been used to argue for 
the presence of planetary mass objects embedded within 
them based on the torque balance gap-opening criterion. 
In the context of this work, we can distinguish four cases: 

i) observed gap with a planet 

ii) observed gap with no planet 

hi) no observed gap and no planet 

iv) no observed gap but with a planet 

The first case is the major topic addressed in this pa¬ 
per, where we suggest a revision of the conditions needed 
to open a gap. Since it may be harder for planets to 
open a gap when considering migration this suggests that 
higher mass planets may be needed to open a gap in 
a disk than that based on the pressure stability crite¬ 
rion of iLin fc Papaloizoul (II993D or the torque balance 
criterion of ICrida et al.l (j2006D alone. Thus the torque 
balance criterion provides a minimum mass for a single 
planet to open the gap. Therefore planets in transition 
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disks may be easier to detect due to the higher expected 
planet masses. In this context it might be worth revis¬ 
iting published results that searched for planets embed¬ 
ded within a gap, or investigating those disks whose gaps 
previously indicated planet masses below the current ob¬ 
servable limits. If this is ruled out, one would be tempted 
then to hypothesize multiple planets. 

The second case would require a more sophisticated 
examination of possible hydrodynamic effects that may 
cause gaps (that is beyond the scope of this current pa¬ 
per). 

The third case, may be common and probably the case 
for the majority of observed T Tauri systems for most of 
their evolution: fewer than 20 % of sun-like stars are ex¬ 
pected to form gas giant pla nets between 0.3 - 20 AU 
based on extrapolations from ICumming et al.l (|2008ll . It 
is important to remember that not all disks are easy to 
resolve in scattered light as observations are limited by 
the obtainable surface brightness. However, if a disk is 
resolved in scattered light and the disk structure does not 
indicate a gap, it is possible that there is no large gap 
of radial extent in the gas surface density distribution as 
the scattered light traces small dust particles which are 
coupled to the gas. But it is an open question whether 
even small amounts of dust might provide enough scat¬ 
tered light to conceal a partially cleared narrow gap in 
the gas. 

The final case, may be more common than current esti¬ 
mates. With state-of-the-art instrumentation, it can be 
difficult to detect a gap that is not completely cleared 
of material unless it is a significant fraction of the cur¬ 
rent beam width of mm-wave interferometers (about 10 
AU). A massive planet with a Hill radius of 5 AU that 
causes a drop in the surface mass density of > 20% might 
just be detectable with high signal-to-noise. Just because 
we cannot detect a gap with current facilitites, does not 
mean that one is not there! In the coming years, ob¬ 
servations with the expanded capabilities of ALMA will 
probe with the best sensitivity and resolution possible 
the presence or absence of gaps in disks around T Tauri 
stars (e.g. the extraordinary images of HL Tau from 
ALMAF^I . Ultimately the next generation of ELTs will 
provide additional constraints on the presence/absence 
of gaps at even higher spatial resolution. 

Finally, since gravitationally unstable disks require 
very high companion masses to open a gap, our results 
suggest that we should not see many gaps (formed by 
planets) in young gravitationally unstable disks. On the 
other hand if gaps in such disks are observed this sug¬ 
gests that some mechanism must take place that slows 
down or stalls planet migration in young disks. Future 
spatially-resolved observations of these disks, e.g. with 
ALMA, will inform us about this. 

6 . CONCLUSION 

We perform 2D hydrodynamical simulations using the 
grid-based code, FARGO, to explore gap-opening in a va¬ 
riety of protoplanetary disks that would be expected to 
form giant planets either by the core accretion or gravi¬ 
tational instability methods. 

The main results of our study are: 

http://www.eso.org/public/news/esol436/ 


• In general, we hnd that gap-opening may be harder 
than t hat predicted by t he torque balance criterion 
alone (iCrida et al.ll2006ll . which is the criterion that 
is commonly used to determine whether gap forma¬ 
tion occurs. 

Typically, this has been considered in the con¬ 
text of stationary planets - however, we find that 
when considering migrating planets, an additional 
timescale factor is very important to consider. We 
find that not only is it crucial whether a planet is 
massive enough to open a gap, but also whether 
the planet can do so quickly enough. We find 
that even if the torque balance criterion is satis¬ 
fied, if the timescale for migration across a region 
in which a planet is opening a gap, tcross, is faster 
than the gap-opening timescale, tgap, a planet is 
not able to open a gap, i.e. for gap-opening we 
require the additional condition Uross > ^gap- Al¬ 
though the timescale asp ect for gap-opening has 
been investigated befor e (|Hourigan fc Wardlll984 
Hdn fc PapaloizoiJfTh^ . we find that the timescale 
comparison might pose a strict limitation to the 
occurence of gap-opening. This is particularly the 
case in migration scenarios where we do not ex¬ 
pect conditions for a ’’smooth” Type I migration 
regime. If a planet’s mass is in the intermediate 
stage (when it is opening a partial gap), it may 
undergo rapid Type Hl-like migration or even run¬ 
away migration. While the torque balance criterion 
may suggest that gap-opening should occur, in real¬ 
ity gap-opening may be prevented by the timescale 
restriction. We show this to be the case with low 
mass disks in which planets are likely to form by 
core accretion, as well as high mass disks in which 
giant planets or brown dwarfs may form via gravi¬ 
tational instability. 

• The disk mass is an important factor that af¬ 
fects gap-opening as this affects the torques on the 
planet and hence its ability to migrate through a 
disk. While the disk mass is not considered in the 
torque balance criterion for gap-opening, it does 
affect the crossing timescale and is thus an impor¬ 
tant factor in the timescale condition, Uross > ^gap, 
which considers all the variables in the torque bal¬ 
ance criterion plus more. 

Our results are applicable to both observational inter¬ 
pretations and theoretical modelling. Since gaps may 
be harder to open for migrating planets than when only 
considering the balance of torques at a particular loca¬ 
tion, the interpretation on companion masses that open 
gaps of particular widths in transition disk observations 
should only be considered as a minimum mass. Further¬ 
more, population synthesis models may find that planet 
survival is more difficult with the addition of a timescale 
constraint. Finally, since the gravitationally unstable 
disk simulations show that planets in such disks seem 
to undergo very rapid migrations, survivability in these 
disks may only be possible if the SG disks are not very 
heavy. Future spatially-resolved observations, such as 
with ALMA, may constrain if gap-opening in such disks 
is possible. 
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APPENDIX 


FURTHER SIMULATION RESULTS 

For the curious reader we present further visualizations of the simulations conducted. Figure lA.ll shows the migrations 
of a 15 Mjup and 45 Mj^p companions in the reference and light SG disks analogous to Figure [21 As shown in Figure 
lA.ll an increased planet mass or decreased disk mass improves the ability for gap-opening. At some stage during the 
inwards migration a gap should open in these simulations (see Figured]) but the very rapid migration in the reference 
SG disks prevents this. In the lighter SG disks the migrations are generally much slower and gap formation appears 
possible. _ 

Figure [AT2] illustrates the evolution of the surface density around an embedded 15 Mjup and 45 Mjup companion in 
the reference SG disk analogous to Figured) The timescale to open a gap (8 - 12 orbits at 70 AU and 20 - 65 orbits 
at 140 AU, w here orbits is defined as an orbital period at 100 AU) is much longer than the crossing timescale (< few 
orbits; Figure |A t] left). As with the results presented in Section 14.1.11 this may be the reason why gaps do not open 
in these disks . 

Figure IA.3I illustrates the evolution of the surface density around an embedded 0.5 Mjup, 1 Mjup and 2 Mjup 
companion at 5 AU and 10 AU in the heavier MMSN-like disk. In this disk the gap-opening timescales are much 
longer than the respective migration timescales (see Fig. [5] on the right), thus preventing gap-opening. 
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Time [Orbits at 100 AU] 


Fig. A.l.— Evolution of the orbital separation of a migrating 15 Mjup (top) and 45 Mjup (bottom) companion in the reference (left) and 
lighter (right) SG disks. The simulation is performed four times each time starting at the same orbital radius but at a different azimuth 
angle leading to four separate migration trails (displayed with different colors). In the reference SG disk the companions migrate rapidly 
through the disk preventing any potential gap-opening. In the lighter SG disk, they decelerate enough or even remain at the radial location 
of introduction so that a gap might evolve. Only the 15 Mj^p reaches the inner grid boundary once, likely due to a stochastic effect. 


Figure IA.4I illustrates the evolution of the surface density around an embedded 30 Mjup companion in the lighter 
SG disk and a 2 Mjup companion in the lighter MMSN-like disk a t 70 AU, 140 AU, 5 AU and 10 AU, respectively. 
A comparison of the gap evolution with Figures [4] (top panel) and IA.3I (bottom panels) shows that the gap-opening 
timescale appears to be independent of mass. 
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Fig. a.2. — Evolution of a gap clearing by a 15 Mjup (top) and 45 Mjup (bottom) companion held on a fixed orbital radius at 70 AU 
(left) and 180 AU (right) in the reference SG disk. For the embedded 15 Mj^p (45 Mjyp) companion a gap forms after approximately 12 
(8) orbits at 70 AU a nd after 65 (20) orbits at 140 AU. The time it takes for a gap to form is larger than the respective migration timescale 
shown in Figure lA.ll for both companions (also see Table [4}. 
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Fig. a. 3.— Evolution of a gap clearing due to a 0.5 Mjup (top), 1 Mjup (middle) and 2 Mjup companion (bottom) held on a fixed 
orbital radius at 5 AU (left) and 10 AU (right) in the heavier MMSN-like disk. For the embedded 2 Mjup companion a gap forms after 
approximately 60 orbits at 5 AU and after 400 orbits at 10 AU. At 5 AU it takes the 0.5 Mjyp (1 Mjup) companion approximately 400 
(200) orbits to form a gap. At 10 AU the 0.5 and 1 Afjup companions do not manage to clear the gap entirely in the 600 orbits simulated. 
In general, the gap-opening timescale for all companions is much larger than the respective migration timescale shown in Figure |5] right 
panel (also see Table Hi. 
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Fig. a.4. — Evolution of a gap clearing by a 30 Mjup companion in the lighter SG disk (top) and a 2 Mjup companion in the lighter 
MMSN-like disk (bottom) held on a fixed orbital radius at 70 AU (top left), 140 AU (top right), 5 AU (bottom left) and 10 AU (bottom 
right), respectively. The gap-opening timescale appears to be independent of disk mass as can be seen when comparing the density 
evolutions here with those in the reference SG disk (Fig [4] top panels) and the heavier MMSN-like disk (Fig. IA.3I bottom panels). 






